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ABSTRACT 

We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space 
Observatory and measure a deficit of surface brightness within their central region after removing 
detected sources. We simulate the effects of instrumental sensitivity and resolution, the source pop- 
ulation, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. 
The amplitude of the central deficit is a strong function of the surface density and flux distribution 
of the background sources. We find that for the current best fitting faint end number counts, and 
excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the 
cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total 
intensity of the CIB of /250/jm > 0.69^Q Q3(stat.)^Q'Jg(sys.) MJy sr _1 , with more CIB possible from 
both low-redshift sources and from sources within the target clusters. It should be possible to observe 
this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which 
would allow tests of models of the faint source component of the CIB. 
Subject headings: cosmic background radiation 

1. INTRODUCTION 

The effect of gravitational lensing is to redistribute the 
intensity from sources behind the lens into images with 
different positions and brightnesses, wh ile conserving 
surface brightness ( Schneid er et al.l I1992T) - This means 
that gravitational lensing not only magnifies the back- 
ground sources, but also changes their apparent density 
on the sky. The details of whether the number counts 
of background sources seen through a foreground grav- 
itational lens such as a galaxy cluster are increased or 
decreased depend on the properties of the lens and the 
slope of the faint counts ()Refregier fc Loeblll997l) , a pro- 
cess known as magnification bias (|Turnerlll980ir 

In this letter, we report detection of a deficit in the sur- 
face brightness of the cosmic infrared background (CIB) 
in the centers of massiv e galaxy clusters m easured using 
the SPIRE instrument (IGriffin et al.ll20Toh o n the Her- 
schel Space Observatory ()Pilbratt et al.l l2010). To inter- 
pret these observations and understand the consequences 
of lensing a background field, the lensing properties of the 
cluster and the background source population are care- 
fully simulated in a large number of realizations. We 
concentrate on the intensity profile after removing de- 
tected sources to highlight faint fluctuations in the CIB. 
The resulting effect after such a removal is a localized 
region of decreased surface brightness at the cluster cen- 
ter. The fluctuations are modulated in our maps due to 
the resolution of the SPIRE instrument, ~ 20 arcsec. We 
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use these observations to constrain the intensity of the 
sub-mm background at 250 /zm. 

2. OBSERVATION OF THE DEFICIT 

This study uses confusion-limited maps of galaxy clus- 
ters from the Herschel Multi-tiered Extragalacti c Sur- 
vey (HerMES: iHerMES Collaboration et aftl20ll . The 
SPIR E data are red uced using a combination of the 
HIPE (|Ott et al.l I2006Q and SMAP (jLevenson et al.ll201d : 
iViero et al.ll2012D packagefl The sample of clusters is 
listed in [HerMES Collaboration et al.l (|2012h and Zem- 
cov et al. (in preparation). 

In order to study the strongly lensed regions it is nec- 
essary to restrict our attention to those clusters large 
enough that f2 c > 3il s , where f2 c is the solid angle of 
the negative magnification region and Q s is the solid an- 
gle of the SPIRE beam. This cut leaves only the four 
clusters from the HerMES sample of twelve: Abell 370, 
Abell 1689, Abell 2219 and RX J 1347-1145. The nega- 
tive magnification region in the image plane maps to the 
interior of the lens caustic in the source plane, so it is 
a reasonable tracer of the area which undergoes strong 
magnification and large deflections by the lensing system. 
The mapping from total mass to lensing characteristics 
is complex and the area fi c is not strictly proportional 
to cluster mass; however, this cut does effectively restrict 
our analysis to the most massive centrally-condensed sys- 
tems. 

To measure the deficit, a catalog is generated for each 
cluster at each SPIRE way elength using the SCAT al- 
gorithm (jSmith et al.l I2012D based on 250 /im source- 
selection. Our goal is not to identify particular point 
sources, but rather to make a catalog which aids in re- 
moving emission from point sources to reveal the dif- 
fuse structure. Thus, employing a low detection thresh- 
old is reasonable. We chose a la threshold, where a is 
the map root-mean-square variation, dominated in these 
deep maps by confusion noise. The measurement is not 
sensitive to the precise value of this number provided it 
is not so high as to leave residual sources (> 3(7 ) nor 
so low as to needlessly mask a large fraction of the map 
(< 0.5tr). This yields a catalog of ~ 500 sources in each 
cluster field which are then subtracted from the cluster 
image using a Gaussian model of the beam. An exam- 
ple of the resulting image for one cluster is shown in 
Fig. [I] which shows a statistical deficit in surface bright- 
ness near the center of the cluster. A simple way to 
visualize this deficit is to plot annular averages of the 
source-subtracted maps, as shown in Fig [5] Similar av- 
erages centered on randomly chosen positions in SPIRE 
images do not show this deficit. 

3. THEORETICAL GROUNDS 

Gravitational lensing conserves surface brightness and 
therefore does not alter the mean intensity of the extra- 
galactic background light at any wavelength. However, 
this statement applies only when all sources are kept in 
calculating the intensity profile and when averaged over 
many directions through the cluster - for some config- 
urations of the background sky unusually bright images 
can be found. The measurement reported here involves 

1 The maps used in this work will be released from the Herschel 
Database in Marseilles (HeDaM, http://hedam.oamp.fr). 




Figure 1. Lensed and field images of the CIB from SPIRE. The 
upper left-hand panel shows the 250 /im Abell 1689 image, and 
the upper right-hand panel shows the same image after detected 
sources have been subtracted. Because the detection threshold is 
equal to the confusion limit and the noise is not uniform over the 
image, at the high contrast level of the right hand image ~ lo" C onf 
noise fluctuations are visible towards the edges of the map. The 
colored contours show the iso-magnification levels from the gravi- 
tational lensing model for a source at z = 1.5 with /i < enclosed 
by the central red contours and fj, = {3, 2, 1.5} in lightening shades 
of blue. The lower two panels show SPIRE 3-band false color sim- 
ulations of the CIB to S < 100 njy, not including noise. The lower 
left-hand panel is a image of the CIB without lensing, while the 
right-hand panel shows the same background that has been prop- 
agated through a lens model for Abell 1689; no emission from the 
cluster or SZ effect is included. 

the removal of bright detected sources, which leads to a 
situation where lensing shows a profile that does not con- 
serve intensity. This does not imply that the total num- 
ber of photons and therefore intensity i s not conserved, 
as discussed in iRefregier fc Loebl (|1997t ). since detected 
sources are removed. 

When viewed through a lensing cluster imposing a 
magnification [i on background sources, a background 
source with flux density S is magnified such that the ob- 
served flux density is S' — The resulting increase 
in the flux density is accompanied by a corresponding 
decrease in the projected surface density of galaxies; the 
observed number density of sources through the lens is 
modified to be N' = iV//i. The overall effect of this 
is to modulate existing fluctuations in the background 
source field. Using the parametrization that the intrin- 
sic number count of the faint background sources scales 
as dN/dS cx S~ T , the differential number counts lensed 
through the cluster become dN'/dS' oc ^ r - 2 S~ r . The 
differential counts imaged through the cluster are cither 
decreased or increased depending on whether the intrin- 
sic counts have a slope T smaller than or greater than 2. 
If sources are not removed fro m the image, the to t al CIB 
surface brightness is conserved IRefregier fc Loebl (|1997lh 

The observations reported here do not constrain the 
exact number counts of sources or the difference of the 
counts of sources through and away from the cluster. 
One reason for this is that the observations are sampled 
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Figure 2. Intensity profiles towards four clusters, illustrating gravitational lensing of the CIB. The left hand plots show the mean flux 
density in 0.25'-wide annular bins for each cluster after all detected sources have been removed. The central data point is constructed 
from all map pixels within the effective radius characteristic of the clusters' critical lines (marked by the vertical dashed lines). At larger 
radii, the data points represent uncertainty-weighted averages starting at the characteristic critical line radius. The mean level of each 
simulated map is constructed to be zero before source subtraction and is reduced by ~ 0.3 MJy sr — 1 after a large fraction of the CIB has 
been detected as sources and removed; the offset varies between targets depending on how many sources are removed. The uncertainties on 
the data points are photometric accuracy of the measurements rather than the precision with which the mean should differ from the model. 
Dotted lines and grey contours show the mean and standard deviation of our sky model, calculated from simulations of the interaction of 
the CIB and each lensing cluster, passed though our data analysis pipeline. The right hand plots show the annular averages of the raw 
cluster images with no source extraction, which highlight the structure of the surface brightness which is subtracted by removing detected 
sources. 



by the limited angular resolution of the SPIRE instru- 
ment, which results in a blending of the faintest sources. 
The effects of limited-resolution observations of the in- 
tensity profile o f background sou rces through the cluster 
are described in lBlainl (| 199 71) and lBlainl ()2002f l. High res- 
olution observations with ALMA will allow one to more 
precisely measure the difference in the counts away from 
and through the cluster. 

4. NUMERICAL SIMULATION 

We can use detailed lensing models of these well- 
studied clusters and models of the source counts and red- 
shift distributions of the sources which comprise the CIB 
to generate numerical simulations of the observed sky 
brightness, and explore how the amplitude of an observed 
deficit depends on source model parameters. To simulate 
the sub-mm backgroun d, we use the model described in 
iBethermin et al.l (|201 If ) . which is tuned to match a vari- 
ety of observed number counts dN/dS a nd is forced to 
integ rate to the total intensity of the CIB (La gache et al.l 
2000). The model associates both a spectral shape and 
redshift to each source, down to z = 0. These simulations 
do not include coherent clustering of the galaxies which 
comprise the sub-mm background, which is a small ef- 
fect on the scale s of interest here (9 < 3 arcminutes; 
iViero et al.ll2012h . Poisson noise is modeled by the sim- 
ulation, and naturally accounts for the cosmic variance 
one would expect in the real sky. The simulated im- 
ages have area 0.25 deg 2 and contain ~5x 10 4 sources 
to source fluxes S m { n = 50 nJy. In the parametrization 
dN/dS oc 5~ r , the faint end number counts of this model 
have r ~ 1.5, so they integrate to a finite value of the 
CIB. 

These simulated backgrounds are t hen lensed using 
lenstool mo dels (jKneib et al.|[l99l I Julio et all 120071: 
UuUo fc Kneibll2009D . which are built using optical mea- 
surements of strong lensing in the cluster fields, mod- 
eling all of the known high-magnification images mea- 



sured for each cluster. The lenstool models replicate 
all of the phenomenology arising from the complex dark- 
matter potentials of each cluster, including multiple im- 
ages, giant arcs, and so on. Gravitational lensing is 
achromatic, so the same lens model applies at all wave- 
lengths. Sources with z < z c are included in the simu- 
lated sky map, but are not lensed; all the sources with 
z > z c are lensed from redshift planes discretized into 
Sz = 0.1 steps. An example of one such simulation is 
shown in the bottom panels of Fig. [TJ Excluding other 
sources of emission, the central surface brightness deficit 
is evident in the modeled image in all three SPIRE bands. 

These simulations show that three factors contribute to 
the effect reported here: (1) the sizes and typical lensing 
amplification factors of clusters; (2) the number densi- 
ties and redshifts of the sources responsible for the CIB; 
and (3) the sensitivity and beam size of SPIRE. These 
combine so that when we examine the lensed cores of 
SPIRE images of cluster fields and remove all reasonably 
significant sources, we are left with deficits in the surface 
brightness of the sky in small regions near the cluster 
center. 

In order to compare the data to the simulations, we 
use the same procedure discussed above on simulations 
of the cluster targets. In these simulations, instrumental 
noise with the same amplitude as measured in the data 
is inserted into the simulated, lensed cluster images, and 
these images are propagated through an identical detec- 
tion/subtraction procedure as the real data, resulting in 
the grey bands in Fig. [2] which represent the mean and 
standard deviation of the model results drawn from 100 
realizations of the background sky that capture the vari- 
ance of the lensing deficit due to different background 
configurations and reflect the range of possible deficit 
shapes and amplitudes arising from our ignorance of the 
exact configuration of the background sources which give 
rise to the deficit. 
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A potential concern is whether simulating sources to 
S ~ 50 nJy is sufficient for the modeling result to con- 
verge. To check this, we have performed simulations 
where sources are not included below S m - ln . By 5 m i n = 
10 /iJy the map is densely populated and the effect has 
converged, and including fainter sources does not have an 
appreciable effect on the amplitude of the deficit. This 
is because sources below 10 /iJy are rarely boosted above 
the la detection threshold and so do not reate a net 
deficit after source removal. 

To investigate the dependence of the deficit on 
the background number counts, we perform simu- 
lat ions drawn f rom the measured counts presented 
m jGl enn et ak| (120101) similar to the calculations in 
iFord et al.l (|2012f >. The slope of the faint end for 0.1 < 
S < 2mJy are varied in the range 1.0 < T < 3.0, 
bracketing the nominal value of F = 1. 65 which ac- 
counts for the full FIRA S background (|Fixsen et al.1 
[l998t ILagache et all 12000). As can be seen in Fig. [3J 
for shallow slopes, the surface density of sources is small 
enough to produce zero flux in the center of clusters. For 
steeper slopes, the probability that a source falls behind 
the center of the cluster image is no longer small and the 
surface brightness increases. Existing constraints on the 
brightness of the CIB at 250 /zm exclude V > 2.3 at 2a 
in this range ([Fixsen et al| [l998). 




Figure 3. Effect of the slope of the faint end source counts V on 
the lensing deficit. The surface brightness of the sky is averaged 
into annular bins centered on the cluster and expressed as a ratio of 
the brightness at radius r to the overall sky surface brightness. For 
slopes similar to that inferred from other measurements (T as 1.5), 
the center of our modeled clusters has an expectation value of zero; 
as the slope is increased, the probability of a source falling into the 
caustic region is increased until the deficit begins to be filled in. 



5. MEASUREMENT OF THE CIB AT 250 MICRONS 

According t o ou r simulations using the 
iBethermin et all ()2011l ) model, statistically the central 
region of these clusters' images have zero surface bright- 
ness from lensed sources. At 250 /im, emission from 
sources within the central arcminute of clusters tends to 
be faint. Therefore a situation arises where the cluster 
itself is invisible, the probability of a foreground source 
coinciding with the cluster center is small, and the 
background is lensed by the cluster, so that no sources of 
emission exist in the center of the cluster. If this region 
has zero surface brightness, as this model suggests, we 
can use it to tie the image to an absolute zero point. 

Of course, if the above assumption is incorrect, any 
sources of emission local to the cluster center - for exam- 
ple, from the BCG or the SZ effect - will bias the mea- 



surement to lower inferred surface brightness. Though 
only a < 0.1% fractional contamination in the 250 /im 
SPIRE band, the SZ effect is significant at longer sub-m m 
wavelengths (jZemcov et al.l l2007: Zc mcov et al.ll2010h so 
we defer reporting the CIB inferred from the two longer 
wavelength SPIRE bands to future work where we model 
the intra-cluster medium (ICM) in each cluster. Sub- 
mm emission associ ated with the centr al region of target 
clusters is known (Ra wle et all I2012T) : for example, in 
the HerMES sample (but excluded from this work) Abell 
1835 and Abell 2390 are both known to host Brightest 
Cluster Galaxies (BCGs) with sub-mm fluxes of many 
mjy. Because of this potential bias, the lensing deficit 
produces a lower limit to the CIB. To check for con- 
tamination from galaxies within the critical region, we 
use Spitzer-MIPS 24 /im data. No BCGs are obvious in 
the SPIRE images of the four clusters, but we can use 
the 24 urn - 250 um BCG flux density ratio measured by 
I Rawle et all ((2012) to estimate an additional one-sided, 
positive-going uncertainty to the central zero point in 
each image. This corresponds to a la value of 0.4 mjy 
per beam at the position of the BCG. We do not detect 
other 24 /im sources associated with the cluster in the 
central deficit region. Diffuse dust emission associated 
with the ICM has never been detected, but from predic- 
tions we expect this kind of emission to be at least an 
order of magnitude smaller tha n the surface brightn ess 
of the lensing deficit at 250 fim (|Popescu et alJl2000T ). 

To determine the absolute intensity of the CIB, the 
average of each map in its central region is used to gen- 
erate a zero point for the image. The mean brightness 
of the resulting image, excluding the central region, is 
then computed. To estimate the uncertainty associated 
with each measurement, we compute the quadrature sum 
of the statistical uncertainty of the maps as traced by 
the standard deviation of the pixels in the central re- 
gion (Stat. 5/250 Aim), the one-sided uncertainty from the 
BCG emission (BCG <5/250/im), the uncertainty from the 
simulations associated with the configuration of the back- 
ground source (Model Wasopm); and the 5 % absolute cal- 
ibration uncertainty of SPIRE (Abs. Cal. 5/250 ^m)- We 
estimate the Eddington bias associated with the effect 
of subtracting sources from the maps before computing 
the mean of the central region by comparing simulations 
including and excluding simulated noise, resulting in an 
estimate of SI < 0.03 MJy sr^ 1 for the Eddington bias on 
the sample. Tablc[T]lists the areas of negative magnifica- 
tion in each target, the inferred CIB brightness for each 
of the four cluster fields that pass all of the selection cuts, 
and the uncertainty budget associated with these mea- 
surements. The statistical-uncertainty weighted mean of 
this value, 0.69^q q? MJy sr _1 , is the inferred brightness 
of the CIB at 250 fim from the lensing deficit method. 
This valu e is consistent with the FIRAS values deter- 
mined bv iPuget et all (| 19961) . iFixsen et all (| 19981) . and 
ILagache et al.l ( 20001 ) in the same band within la. 



6. DISCUSSION 

We have detected an interesting new phenomenon 
which is due to the properties of the sub-mm background 
and gravitational lensing in massive clusters. This deficit 
in the measured surface brightness constrains the small- 
est allowable surface brightness of the CIB and so is an- 
other way to limit its absolute brightness, adding to a 
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Table 1 

Lensing deficit CIB limits & error budget. 



Cluster 




CIB ^250 Mm 


Stat. 5/250 urn 


Model <5/ 2 50 Mm 


BCG 5/250 Mm 


Abs. Cal. 5/250 Mm 


Total 5/250 Mm 




(«.) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


Abcll 370 


3.0 


0.67 


±0.08 


±0.11 


±0.05 


±0.03 


-0.14, +0.15 


Abell 1689 


7.5 


0.68 


±0.06 


±0.09 


±0.05 


±0.03 


-0.11, ±0.12 


Abell 2219 


3.3 


0.71 


±0.09 


±0.08 


±0.05 


±0.04 


-0.13, ±0.14 


RX J 1347-1145 


11.3 


0.71 


±0.05 


±0.09 


±0.05 


±0.04 


-0.11, ±0.12 


Total 


25.1 


0.69 


±0.03 


±0.05 


±0.09 


±0.04 


-0.07, ±0.12 



list which includes P{D) analyses, stacking, and source 
counting, in addition to direct photometric measure- 
ment. Furthermore, in principle the presence of this 
deficit places interesting constraints on the faint end of 
the number counts, although larger samples and compre- 
hensive modeling are required to convert measurements 
to source count constraints. 

Because gravitational lensing is achromatic, this effect 
occurs in all of the SPIRE bands. Other instruments 
working at sub-mm and mm wavelengths like ACT/SPT 
and CCAT should also be able to measure this deficit 
effect, assuming that the background in question is suf- 
ficiently far behind the lensing cluster, and that the in- 
strumentation is sensitive to sources at levels similar to 
SPIRE's with enough angular resolution to resolve the 
central region of clusters. In addition, this effect pro- 
duces a complicated, spatially structured CIB surface 
brightness distribution which is a potential foreground 
for high resolution SZ effect measurements at longer 
wavelengths. 
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